Li electron spectrum versus neutrino spectrum: implications 
for the Sudbury Neutrino Observatory 



G. Jonkmans^'^*, I.S. Towner^ and B. Sur^ 

Department of Physics, Queen's University, Kingston, Ontario, Canada K7L 3N6 
^ AECL, Chalk River Laboratories, Chalk River, Ontario, Canada KOJ IJO 

(February 9, 2008) 



Abstract 



The sensitivity of the Sudbury Neutrino Observatory (SNO) to measure 
the shape of the recoil electron spectrum in the charged-current reaction of '^B 
solar neutrinos interacting with deuterium can be improved if the results of a 
^Li beta-decay calibration experiment are included in the test. We calculate 
an improvement in sensitivity, under certain idealistic assumptions, of about 
a factor of 2, sufficient to resolve different neutrino-oscillation solutions to the 
solar-neutrino problem. We further examine the role of recoil and radiative 
corrections on both the neutrino spectrum and the ^Li electron spectrum 
and conclude that the influence of these effects on the ratio of the two spectra 
as measured by SNO is very small. 
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I. INTRODUCTION 



The Sudbury Neutrino Observatory (SNO) |1| will utilize the interaction of *B solar 
neutrinos with deuterium in heavy water to measure the shape of the recoil electron spectrum 
in charged- current (CC) interactions and the ratio of the number of charged-current to 
neutral-current (NC) events. In what follows, we will concentrate only on the test of the CC 
spectral shape, where one of the measurable quantities is the average kinetic energy of the 
detected recoil electrons, (Te)y, to be defined more precisely below. The main uncertainties, 
other than counting statistics, arise from uncertainties in (a) the theoretical standard 
neutrino spectrum, (b) the detector energy resolution and (c) the detector absolute energy 
scale. If the uncertainties from these sources can be reduced sufficiently SNO should be 
able to distinguish between a Standard Solar Model with no neutrino oscillations and one 
with various choices of neutrino oscillation scenarios [0: a Small Mixing Angle (SMA), a 
Large Mixing Angle (LMA) or Vacuum Oscillations (VAC) solutions. One way to reduce 
these uncertainties is to introduce a *Li calibration source. As part of its overall calibration 
strategy, the SNO collaboration will install a system that is capable of producing a ®Li 
source placed at several different locations inside the detector 0. Detection of electrons 
from this source would demonstrate that the results obtained by SNO for a known beta- 
decay spectrum are consistent with those measured in the laboratory. 

In this paper we point out that uncertainties in the measured CC spectral shape can 
be reduced considerably by a direct comparison with the measured beta-decay spectrum 
of ^Li. For instance, rather than considering {Te)y to be the benchmark for comparing 
experiment with different theoretical expectations, we consider the ratio {Te)u/{Te)e) where 
(Te)e is the average kinetic energy of detected electrons from ^Li decay. In forming ratios 
such as this, many uncertainties cancel producing a much improved error budget. This 
suggested analysis strategy is based upon the fact that the electron spectra from ®Li beta 
decay and from solar-neutrino absorption on deuterium originate from mirror weak decays 
to common final states in *Be*, hence are highly correlated, have somewhat similar shapes 
and cover essentially the same energy range. In formulating this strategy, we have calculated 
the forbidden and radiative corrections to the '^Li beta and ^B neutrino spectra and have 
considered the effects of theoretical and experimental uncertainties on the above-mentioned 
ratio. 

In Sec. H we define more carefully the quantities involved and outline the suggested 
strategy. In Sec. |T| we give some calculated results, and in Sec. ^ we summarize our 
findings. Most of the formulae have been relegated to appendices. 



II. THE MEASURED QUANTITIES 

The recoil electrons produced by ®B solar neutrinos being absorbed on deuterium have 
an average kinetic energy (for an ideal detector) given by 

,y .ideal ^ /t,.„ dT^T^ J dEA.{E.)PUE.r-^{E^) 

where T^[^ is the threshold kinetic energy below which events in SNO will be discarded. For 
the discussions that follow, we adopt Tmin = 5 MeV. The charge- current differential cross- 
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section for neutrino absorption on deuterium dacc/dTg, is a function of neutrino energy. 
We will use the calculated cross-section of Ellis and Bahcall [Q for our computations, and 
adopt from Bahcall and Lisi a 0.14 % la error in {Te)y, this being the difference when 
Kubodera-Nozawa cross-sections are used instead. The neutrino spectrum is denoted 
\u{Ey) = dT/dE,^, and Pee{Ei,) is the probability that an electron-neutrino produced at 
the core of the sun remains an electron-neutrino by the time it is detected on earth. The 
different neutrino-oscillation scenarios yield different survival functions, Pee{Ev)- In the 
Standard Model with no neutrino oscillations, Pee{Ey) = 1. 

The decay electrons produced in ^Li beta decay have an average kinetic energy (for an 
ideal detector) given by 



i ideal It^;„ d-TeTeXe(Te 



where Ae(T'e) = dT/dE^ is the electron spectrum. The beta decay of ^Li is the isospin 
analogue of the beta decay of ^B; both populate the same daughter ^Be states. The ^Be 
excited states are unstable and break up into two alpha particles. As a consequence, the 
shape of the ^Li (3~ spectrum and of the ^B neutrino spectrum deviate significantly from the 
standard allowed shape. Measurements of the delayed a spectra allow one to determine 
the profile of the intermediate ^Be state and thus to calculate the deviations in the electron 
and neutrino spectra. This is discussed at length in Bahcall et al. |P, where it is shown 
there is considerable discrepancy among different experiments related to the absolute energy 
calibration of the alpha particles. Bahcall et al. choose to display this uncertainty as a 
possible offset h in the energy of the alpha particles from ^Be break up: E^ —>■ E^ + b. The 
value of the offset affects the calculated shape of both the beta and neutrino spectra. The 
effective 3a uncertainty of the offset is estimated to be 6 = ±0.104 MeV this estimate 
includes theoretical errors. Bahcall et al. have tabulated their recommended ^B neutrino 
spectrum together with their 3a errors. 

In this work we require both a ^B neutrino spectrum and a ^Li electron spectrum from 
decay to the same final states in ^Be*, in order to compare the effects of common-mode 
uncertainties. Therefore we use R-matrix theory to model the profile of the intermediate 
^Be state, and in particular, the R-matrix fits of Barker . To reproduce the possible offset 
in the energy of alpha particles, we shift the energy of the dominant resonance, Ei — >■ Ei + b', 
where the offset b' is obtained from the requirement that this method gives the same error as 
the Bahcall et al. treatment. The same offset is used consistently for both decay spectra. 
Details of the R-matrix analysis are given in Appendix 0. 

There are two further differences between the ^B neutrino spectrum and the ^Li electron 
spectrum originating in (a) recoil (or forbidden) corrections and (b) radiative corrections. 
Details of these have been relegated to Appendices P and |C| respectively. However, in Figs. |l| 
and we show plots of these corrections for the individual spectra and their ratio. Although 
both corrections are energy- dependent and hence a factor in determining (Tg) their ratio is 
very much less energy dependent. Thus any error in (Te) associated with determining these 
corrections is much reduced in considering the ratio (Tg)^/ {Te)e. 

In Fig. |I] we show the recoil correction for the ®B neutrino spectrum and ®Li electron 
spectrum for two cases. The calculation is based on the elementary-particle treatment of 



beta decay of Holstein ||T0[ and the correction depends principally on two parameters: b/c 
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and d/c. Here 6, c and d are the weak-magnetism, Gamow- Teller and induced pseudotensor 
current form factors respectively. The two cases correspond to the induced pseudotensor 
current being retained in or removed from the calculation. In both cases the ratio is almost 
independent of energy, indicating little sensitivity to the presence of induced pseudotensor 
currents. 

In Fig. 1^ we show the radiative correction for the neutrino spectrum and ^Li electron 
spectrum again for two cases. In the first case, we assume the internal bremsstrahlung 
photons are not detected and so the energy of the photon is integrated over in obtaining 
the radiative correction. In this case there is a large difference in the energy dependence 
between the radiative corrections for the electron and neutrino spectra. In the second 
case, we account for the fact that in a calorimetric detector such as SNO, the internal 
bremsstrahlung photons from a ®Li source placed in the detector would be recorded and 
summed with the beta energy. We assume, for simplicity, that the efficiency and detector 
response for electrons and photons is the same, and compute the radiative correction for the 
summed energy (photons plus electrons) deposited. Now one sees there is a much reduced 
difference between the energy dependence of the radiative corrections for the electron and 
neutrino spectra. 

From Figs. |I| and|^, it is fairly clear there will be little uncertainty in the ratio, {Tf.)^/ (Te)e, 
arising from any uncertainty in recoil and radiative corrections. Therefore, we will not con- 
sider these corrections explicitly any further; but rather we will assume they are effectively 
absorbed in the energy offset parameter, 6, which as has been mentioned above was estimated 
by Bahcall et al. |§] to be ±104 keV, and contained theoretical uncertainties associated with 
recoil and radiative corrections. 

Finally, we turn to detector-related uncertainties. The measured electron kinetic energy, 
Te, determined by SNO will be distributed around the true energy, T^, with a width es- 
tablished by the detected Cerenkov photon statistics. The resolution function R{T^, Tg) is 
approximated by a normalized Gaussian: 



with (Tio the resolution width at = 10 MeV. We follow Bahcall and Lisi and use 
o"io = 1.1 ±0.11 MeV {la errors) as an illustrative example. 




(3) 



where a = cr(Tg) is an energy- dependent la width given by ^ 




(4) 




(5) 




U^^^.^JT,TJdT^R{Ti,T,)X,{Z) 



(6) 
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The other major uncertainty concerns the absolute energy cahbration for the SNO detector. 
For this, we assume a 1% error at 10 MeV. This can be easily implemented in the calculations 
by modifying the energy- resolution function, Eq. (^, with the replacement 

i?(T,',T,)^i?(T,' + 5,T,), (7) 

where 6 = ±100 keV (la error). 

The above methods for treating detector-related uncertainties are a rough parametriza- 
tion of detailed considerations regarding the spatial and directional response to energy de- 
position in the SNO detector. Other sources of uncertainties which have not been taken into 
account in this work are (a) the energy dependence of the detector efficiency, (b) the statisti- 
cal and systematic uncertainties associated with background subtraction and neutral-current 
event separation, (c) the uncertainties in the ^Li beta-spectrum (as measured in SNO) due 
to the source container, and (d) the position dependence of the measured energy. Evaluation 
of the effects of these uncertainties will have to await detailed Monte-Carlo simulations and 
actual measurements of detector performance when SNO is operational. 



III. THE ERROR BUDGET 



Having listed the source of uncertainties, it remains to quantify their contribution to the 
error budget for the detection of recoil electrons following neutrino absorption on deuterium 
alone, (Tg)^, or in concert with a calibration experiment with a ^Li source in the detector, 
(Te)i,/(Te)e. The results are given in Table |. The error budget for the former case has 
been given by Bahcall and Lisi ^ and we follow their example. The statistical error {3a) 
on 5000 CC events is estimated at 0.98% for {Te)u- In the calibration experiment there will 
be considerably more events counted, so the statistical uncertainty in the ratio {Te)u/ {Te)e 
will be dominated by the 5000 CC neutrino events: we use the same statistical error. For 
the neutrino-absorption on deuterium, the uncertainty in the absorption cross-section is 
common to both {Te)u and {Te)u/ {Te)e, so again we use the same error. The uncertainty in 
the neutrino spectrum from the beta decay of is characterised by the offset parameter, b, 
associated with the absolute energy calibration of the a-particles in the measured delayed a- 
spectrum. Upto differences associated with isospin-symmetry breaking, recoil and radiative 
corrections, the same uncertainty occurs in the calibration experiment, so the error in the 
ratio (Te)j,/(Te)e is much reduced. Table | shows a reduction of uncertainty by a factor 
of 7.6. The same reasoning follows for the detector-related uncertainties, since in both 
experiments, energetic electrons are being counted in the same detector. We find a reduction 
in uncertainty by a factor of 8.5 due to the detector's energy resolution, and a factor of 5.0 
due to its absolute energy calibration. The energy calibration of the detector here has been 
oversimplified. The analysis above assumes that the electron sources from ^Li beta decay and 
from neutrino absorption are identical, but of course they are not. The ^Li calibration sources 
are highly localized in space and time, while the recoil electrons produced by absorption of 
solar neutrinos are distributed throughout the detector and distributed in time. So our 
analysis here represents an optimal scenario, but it is sufficient to make the point that 
an improvement in the error budget will be achieved by performing an in situ calibration 
experiment. 
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If all errors are added together in quadrature, then the 3a uncertainty in the CC test 
at SNO using (Tg)^ alone is 2.8% and is dominated by the uncertainty in the absolute 
energy calibration. However, in combination with a ^Li calibration experiment, the 3a un- 
certainty in the ratio (Te)y/(Te)e is reduced to 1.2%, and more importantly, is dominated 
by counting statistics. This analysis strategy will significantly improve SNO's ability to dis- 
criminate among the different neutrino-oscillation scenarios. In Fig. ^ we display the error 
budget together with the theoretical value of {^e)v/ {Te)e for various neutrino-oscillation so- 
lutions to the solar-neutrino problem: the two (best-fit) Mikheyev-Smirnov-Wolfenstein 
(MSW) solutions at small and large mixing angle (SMA and LMA) and the purely vacuum 
(VAC) oscillation solution. In an analagous plot by Bahcall and Lisi |^, it is shown that 
a measurement of {Te)^ alone is unlikely to resolve the two MSW solutions. On the other 
hand, the ratio measurement. Fig. |^, clearly shows the LMA and SMA solutions as being 
experimentally distinguishable, if the other uncertainties in the detector can be minimized. 

Although we have considered the ratio of (Tg)^/ {Te)e, it is possible to construct other 
experimental quantities from the measured CC and ^Li spectra. One of the simplest observ- 
ables to consider is the actual ratio of the normalized measured CC spectrum, A^!^, to the 
normalized measured '^Li electron spectrum, X^^p. Because of the high degree of similarity 
between those two spectra, the effect of uncertainties on the detector response function are 
expected to largely cancel when forming this ratio. In Fig. ^ we display the effect of the 
uncertainty in the detector energy scale, 6 on the the double ratio 



where Xtheor refers to a spectrum where all parameters including the absolute energy calibra- 
tion are assumed to be known. If Af'^^ is not measured (i.e. the denominator is assumed to be 
unity), then the dashed lines enclose the error band with ±3a errors on 6. A measurement of 
Af^p allows the cancellation of correlated errors and produces the much reduced "±3cr" error 
band enclosed by the solid lines. The effect due to the statistics of 5000 CC events above 5 
MeV is shown by overlaying on Fig. ^ the shape distortion expected for the best-fit mass and 
mixing values for the SMA solution (Am^ = 5.4 x 10"^ eV^ and 8111^26 = 7.9 x 10"^). 
We can clearly see that making use of the ®Li spectrum as measured in SNO reduces the 
effects of systematic uncertainties and greatly increases the discriminating power of the CC 
shape measurement. 



The SNO detector is being constructed primarily to measure the neutral-current (NC) 
to charged-current (CC) cross-section ratio for solar neutrino absorption in deuterium. It 
is anticipated this measurement will give a strong and unambiguous signal of neutrino os- 
cillations (or the lack thereof). Another test, which does not rely on the knowledge of the 
absolute reaction cross-sections, is a measurement of the spectral shape of the CC reaction 
alone. To improve the sensitivity for this measurement, a calibration experiment is planned 
at SNO with a *Li source placed in the detector. If instead of considering only the first 
moment of the recoil electron spectrum in the CC reaction, {Te)u, the ratio of this moment 



R{CCrLi) 




'expl '^theor 



(8) 



IV. SUMMARY 
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to the equivalent moment in the ^Li beta decay is used as a benchmark, then an increase in 
the sensitivity of a factor of 2.4 could be achieved in the CC-shape test. This improvement 
would be sufficient to resolve the small-angle and large-angle MSW solutions to the solar- 
neutrino problem. We stress again, however, that not all uncertainties are known at the 
present time and some will require actual measurements when SNO is operational, so our 
analysis here represents an optimistic best-case scenario. To us, a calibration experiment 
with a ^Li source certainly seems beneficial. 
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APPENDIX A: R-MATRIX FIT 

The allowed /5-transitions from the decay of ®Li or ®B populate 2"*" states of ^Be, which 
then decay into two a-particles. The a-spectra show a pronounced peak, corresponding to 
the 2"*" first excited state of ^Be at ~ 3.0 MeV. Attempts to fit the a-spectra assuming 



only one state in Be is operative fail to give enough yield at high energies. Barker |T2 



proposed R-matrix formulae in the many-level, one-channel approximation and used them 



to fit the /5-delayed a-spectra. More recent fits have been given by Warburton |T^ and 
Barker |]. 

In the many-level, one-channel approximation the beta-decay differential cross-section is 
written 

dT oc -g{E,)peE,ElF{±Z, E,)5{Q,, + m - E., - E, - E,)dE,dEJE,, (Al) 

where Qec + m = M — M', with M the mass of the parent nucleus, M' the ground-state 
mass of the daughter nucleus and m the electron mass. Here E^ is the excitation energy 
in the daughter nucleus, E^ the electron energy, the electron momentum and F{Z, E^) 
the usual Fermi function. This expression differs from the standard one by the presence of 
the function g^Ex)/^ and an additional integration, dE^^, over the excitation energy of the 
daughter nucleus. 

The function g{Ex) in R-matrix theory is 

I R,{E,) P P^iE,) 



^^^^^ \l-{S2{E,)~B2 + iP2{E,))R{E,) r 
R2{E.)=Y. 



Jx{Mgt)x 
T E,-E ' 
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where E is the channel energy, i.e. the energy above the a + a threshold: E = Ex + 0.092 
MeV. Here P2{E) and 5*2 are the penetrability and shift factors for L = 2 partial waves 
and are expressed in terms of Coulomb functions for a + a scattering evaluated at a chosen 
channel radius, 02- Finally, B2 is the boundary condition parameter. The sum over A is 
a sum over all resonances retained in the calculation. For each resonance there are three 
parameters: Ex the centroid energy of the resonance, 7a the reduced width amplitude for 
a + a scattering, and {Mgt)x the Gamow- Teller matrix element for the beta-decay feeding 
of the resonance. 

In the limit of a single narrow resonance, 7A-f2 ~* 0, then 

g{Ex) ^ {MGT)l7r6{Ex - E) (A3) 

and the rate expression, Eq. reduces to the standard one. To obtain the electron 

spectrum, we integrate Eq. (^T|) over neutrino energies 

dV 1 rQec+m-E, 



oc peEeE{±Z, Ee)- / dExg{Ex){Qec + m-Ex- EeY, (A4) 

TT Jo 



dE, 

while for the neutrino spectrum, we integrate over electron energies 
dT 1 pQec-Ey 

-—(xEl- dExg{Ex)(3{Qec + m - E, - E,fE{±Z, Qec + m-E,- E,), (A5) 

arji, 71 Jo 

where f3 = pjE, = [1 - mV(Qec + m - E, - E,ff\ 

It remains to specify the parameters, which we take from the work of Barker |^ , where fits 
are made simultaneously to the a+a scattering (i-wave phase shift as well as the beta-delayed 
a-spectra. Separate fits, however, are done for the ®B and ^Li beta-delayed a-spectra. Thus 
there are differences in the parameter sets between the two cases. However, these differences 
are small and can be attributable to isospin-symmetry breaking. The adjustable parameters 
in the formulae include the channel radius 02, the boundary condition parameter i?2, and the 
eigenenergies Ex, reduced width amplitudes 7a and Gamow- Teller matrix elements {Mgt)\ 
for the various 2"*" levels, A. The best fits were obtained with a large channel radius of 
02 — 6.5 fm. For this value of 02, the second 2"*" level is at about 9 MeV, with a width 
of about 10 MeV. To date, there is no evidence for a state in this vicinity. Altogether 
five resonances are included, denoted A = 1,2,3,0,1'. The state, A = 1, is the dominant 
resonance corresponding to the first excited state in ^Be at about 3 MeV; while A = 2 is 
the second 2+ state just mentioned, and A = 3 represents a background state well above 
the energy range being fitted and which naturally is not fed in beta decay, {Mgt);^ = 0. In 
addition, there are two narrow states at Ex = 16.6 and 16.9 MeV, which are isospin mixtures 
of T = and T = 1 2+ states, the latter being the analogue of the ground states of ®Li 
and ®B and is fed by Fermi beta decay, {Mp)i/ = \pl. The method of handling this isospin 
mixing is described in Barker and Warburton |T3| . 



The fitted parameters Ex., 7a and {Mgt)\ depend on the choice -82- Identical fits can 
be found for different choices; and connection formulae are available []T1| to relate one set of 
fitted parameters to another. The standard choice is to set = 5*2 (-E^), where is one of 
the resonances of the fit. The parameters obtained with this choice are labelled E^^\ '-/^^^ and 
{Mgt)[^\ Barker § only gives the parameter values for the case /i = A, and the connection 
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formulae have to be used to relate them all to a common B2 value. This is cumbersome to 
use; for our purposes it is sufficiently accurate to take the parameter values for = A and to 
vary B2 as follows: for < 5 MeV where the A = 1 resonance dominates, use B2 = 6*2 (-Ei); 
for 5 < £"3; < 13 MeV where the A = 2 resonance dominates, use B2 = 6*2 (-E2); and for 
Ex > 13 MeV where the doublet states dominate, use B2 = 5*2 (-E), where E is the average 
energy of the two isospin-mixed states. The parameter values are given in Table |I[ 

APPENDIX B: RECOIL CORRECTIONS 



Recoil corrections to allowed beta decay have been given by Holstein [10| and used by 



Bahcall and Holstein |jT5[ in their discussion of corrections to the spectrum of neutrinos 
produced by the beta decay of in the sun. There are two sources to the recoil corrections: 
(a) true kinematic corrections arising from relaxing the approximation that the recoiling 
nucleus is at rest, and (b) the introduction of induced terms into the V-A weak hadronic 
current, mainly the induced tensor term in the vector current (weak magnetism) and the 
induced pseudotensor term in the axial current. We have not included an induced scalar 
term in the vector current, as it gives a zero contribution under the assumption of the 
conserved vector current (CVC) hypothesis, or an induced pseudoscalar term in the axial 
current as it is small in beta decay. The beta decay spectrum, therefore, is predominantly 
given by four form factors, denoted in the formalism of Holstein |]10| by a, b, c and d for 



the Fermi, weak magnetism, Gamow- Teller and induced pseudotensor current form factors 
respectively. 

Let P, pr, Pe and p^ denote the respective four momenta of the parent nucleus, daughter 
nucleus, electron and neutrino. Further, let M be the mass of the parent nucleus, MAr(^Li) or 
M^i^B); M' be the ground-state mass of the daughter nucleus, M^r C^Be); E^ the excitation 
energy above the ground state; and m the mass of the electron^. Then we define 

q = P -Pr=Pe+Py, 

A = M - M' - E^. (Bl) 

Further the maximum electron energy is 

Ao = (M=^ + - M'^) /2M 
/\? -rri^ 

~ A , B2 

and the maximum neutrino energy is 

Co = [m^ - (m + M'f] /2M, 
C..™.A-<^. (B3) 



^Holstein |1C] writes Mi and M2 as the masses of the parent and daughter nuclei respectively, 



and M as their arithmetic mean. 
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The form factors are functions of the four-momentum transfer squared. It is convenient to 
expand these form factors 



a(g2) = ai + a2iq^/M^) + ... 
c(g2) = ci + C2(gVM2) + ..., 
b{q') = b + ..., 
d{q^) = d + . . . 



(B4) 



and retain just the terms shown. In the formulae that follow, we will drop the dependence on 
the four- momentum transfer, i.e. neglect 02 and C2, and not display the electromagnetic cor- 
rections, which are of order {aZ). Both these effects are very small on the recoil corrections, 
although they are retained in the final computations. 

The probability that an electron of energy. Eg, is produced in an allowed transition in 
beta decay is proportional to 



dT 
dK 



(B5) 



where re{Ee) is the recoil correction, F{±Z,Ee) is the usual Fermi function and pe is the 
electron momentum, pe = [El — rn?Y^'^ . The upper sign is for electron emission in /3~ decay, 
the lower sign for positron emission in (3^ decay. The recoil correction is 



' ^ 2 A{cl + cid± cib) ^2Ee (Saj + bcj ± 2cib) 



3M 



(«i + cf ) 



3M 



1 {2cl + Cid±2cib)){Ao- Ee)' 



3 ME, 



[aj + cj] 



(B6) 



This result has been given before by Holstein ||10 |. 

For the neutrino spectrum, the probability that a neutrino of energy, E,^, is produced in 
an allowed transition in beta decay is proportional to 



dV 

dK 



oc {aj + cl)F{±Z,A- E,)El{A - E,)[{A - E, 



m 



211/2 



(B7) 



where r^{Ei,) is the recoil correction given by 



1 - 



rriA 
Id 



A- E^ A- E^ + m 



2 A (c? + cid T cib) 2 E^ {Saj + hc\ ^ 2ci6) 



+ 



3M 
1 



3M 



{2c\ - cirf T 2ci&) I (Co + m - E,)\{C^ + m - E,f - m^/^ 



[al + cl) 



j (A-E,)[(A-E,)2-m2]V2 



fB8) 



Note that r,j{Ey) is not given by rf.{E(. A — E^). The reason is that, in recoil order, 
the energy available is distributed three ways: to the electron, to the neutrino and to the 
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recoiling nucleus. The replacement ^ A — Ei^ is only correct in the approximation that 
the recoiling nucleus is at rest. 

It remains to specify the values for the form factors. For the dominant transition in the 
decay of or *Li, 2"*", T = 1 — 2+, T = 0, the transition is pure Gamow- Teller and ai = 0. 
Then, for the recoil correction it is only necessary to supply the ratios b/ci and d/ci. For 
these we use the same values as Bahcall and Holstein ||T5| : 

b 

7.7 ± 1.0, 



Aci 

-^ = 1.9 ±1.3, (B9) 

Aci ^ ^ 

where A is the mass number, A = 8. The weak magnetism value is determined from the 
CVC hypothesis and the measurement of the Ml width of the 2+ analogue state in ^Be [|16 



This value is consistent with, but less precise than, the recent measurement of the Ml width 
of de Braeckeler et al. The induced pseudotensor form-factor value is determined from 



fits to the measured j3-a correlations on ®B and Li [|1^. The form factor d is the lesser 
precisely known of the two listed in Eq. (|B9D . To sample the dependence on ci, some results 
will be given for ci = 0. 



APPENDIX C: RADIATIVE CORRECTIONS 

Radiative corrections to the electron spectrum from allowed beta decay have been con- 
sidered in a number of papers [p!9|-23|. In obtaining the corrections, integrations are carried 



out over the allowed neutrino and photon energies and the results exhibited as a differen- 
tial spectrum in electron energy. The contributions to the radiative corrections have two 
components: the emission of real photons (internal bremsstrahlung) and virtual radiative 
corrections due, for example, to the exchange of photons between charged particles. The dif- 
ferential rate for real photon emission in an allowed beta transition is given by the expression 



24 



a 1 

dV (X - — Q^dQEldE^peEedEedx5{A - Ee- E^-e) 
'Ee + e\ /32(l-a;2QVe2) ^ ^ 



Ee J {e-pQxf El{e-l3Qx) 



(CI) 



Here e represents the photon energy, e = [Q"^ + A^]^/^, Q being the photon momentum and 
A a small nonzero photon mass introduced to regulate the infrared divergence. Further, E^. 
is the electron energy, the electron momentum, (3 = Pe/E^, A the maximum electron 
energy and x the cosine of the angle between the electron and photon directions. The delta 
function is used to integrate over the neutrino energies. The integrations over Q and x are 
however very delicate: the logarithmic pole in A has to be extracted to cancel with the A 
dependence coming from the virtual corrections. 

In the calibration experiment at SNO, a '^Li source will be placed inside the SNO detector. 
The real internal bremsstrahlung photons emitted by this source will in principle be detected. 
So we cannot follow the normal procedure of obtaining the radiative correction by integrating 
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over Q. Suppose that SNO detects photons of energy greater than some threshold, say uj. 
Then for Q less than tu, the normal procedure can be followed except that in the integration 
over Q, the upper limit Qmax is taken to be the lesser of or A — i^e, rather than just 
A — E'e- The result for real photon emission is 



-— (X (A - E,yp,E,-g^{E,,A) 
ahip vr 



^<(^e,A) 



Ely 



rtii 



2/3 



In 



y^E, 
+2 In 



1 



A 



2/3 



In 



,1-/5. 



(C2) 



where y is (A — i?e), and Qmax the lesser of u or y. The value of C is given by Kinoshita 
and Sirlin 



C = 21n2 



2/3 



In 



,1-/5. 



where L(2) is a Spence function 



L(z)= r^ln(ll-tl). 
JO t 



,1-/5. 
'1 - 



2 + ln 



/?^ 



L 



(C3) 



(C4) 



For the case when the photon energy exceeds u we must proceed a little differently. We 
assume, for simplicity, that photons and electrons are detected with equal efficiency and 
that it is the sum of the energy deposited that is observed in the SNO detector. Thus in Eq. 
( |U1| ) we change variables from dE^dQ to dEedX, where X is the sum Eg + Q, and integrate 
over Ep from m to X. The result is 



^^iA-Xf/3iX)X'-g^{X,A), 

dX TT 



g>{x,A) 



/3iX)X 



F{X) In 



X — m 



/3iX)X 



^^^FiX)-FiE) 
X-E 



1 (XJi - J2) 

2 /5(X)X2 



(C5) 



where 



■h 
J2 

F{E) = /3E 



\l-/3 

X 

dEE In 



2/3 V 1 



.1-/5. 
l + /?\ 



P, 



(C6) 



12 



and p = f3{E) 



m 



'■ /E'^y/'^. To these expressions must be added the contribution from 



virtual radiative corrections, which have been given by Yokoo et al. 

oc (A - 



The result is 




9v{Ee,A) 

A 



(C7) 

where M is the nucleon mass, and A a renormalization scale. In the early years, Yokoo et al. 



22 invoked an intermediate vector-boson model to argue that A should be of the order of 



the nucleon mass. Following the development of the Weinberg-Salam Standard Model, the 
virtual radiative correction includes additionally the exchange of Z-bosons. Sirlin |^ has 
shown, remarkably, that the form of the expression, Eq. ( |C7|) , remains unchanged except 
that A is now replaced by the mass of the Z-boson, m^. For our purposes, the value of 
A is not important as it only enters in a constant, energy- independent term and would be 
absorbed into the normalization of the spectra. In Fig. ^ we used A equal to the nucleon 
mass, M. Note that the infrared divergence term in ln(A) exactly cancels between the real 
and virtual photon expressions, as it should. In summary, the radiative correction to the 
electron spectrum is 

dV 



^^iA-Eff3iE)E'R,iE,A), 
R,{E, A) = 1 + - (g^iE, A) + g,{E, A)) 

TT 

= l + -{gy{E,A)+g,{E,A)) 



for E 



Ee < UJ, 



for E = X> UJ. 



(C8) 



Finally, we consider the radiative correction for the neutrino spectrum. This is not 
obtained from the electron spectrum radiative correction by the substitution oi Ee ^ A — E^ 
because the energy available to the leptons is distributed three ways: to the electron, to the 
neutrino and to the internal bremsstrahlung photon. This was first pointed out by Batkin 
and Sundaresan . Since the neutrinos are originating from beta decays in the sun, there is 
no question this time of the internal bremsstrahlung photons being detected. Thus, starting 
from Eq. ( |CT| ) for the real photon radiative correction, the delta function is used in the 
integration over electron energies, while the integrations over Q and x are done carefully to 
isolate the infrared singularity. The result has been given by Batkin and Sundaresan . 
The virtual radiative correction is unchanged and given by Eq. Putting these together 

we get 



dV 
dE~, 
RAE..A) 



<^El{A-E,f(3Ru{E.,A), 



l + -{g,{E,,A) 

TT 



g,{A-E,,A)) 



g,iE„A) = 2\n 



A-E„ 



m 




C + 



1 



2(A - E,)^f3 



(C9) 
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where 







A-E^-m Q 

dQ ^ 



2/3(A - E,) 



{A-E,-Qf-m^V^'F{(3m 



2(3{Q) [l-P{Q)) 



m) 

and /3 = /?(0). 



1 - 



{A-E,-QY 



1/2 



(CIO) 
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TABLES 



TABLE L Percentage 3a error in (Tf,),/ from Bahcall and Lisi ||5) and in the ratio, {Te)^/ {Te)e- 



% Error in: {T^ {T^)J{T,)^ 

Statistics of 5000 CC events 0.98 0.98 

CC cross-section 0.43 0.43 

Neutrino spectrum 1.14 0.15 

Energy resolution 0.94 0.11 

Absolute energy calibration 2.04 0.41 

Total error 2.82 1.16'" 



^We stress this is an optimistic estimate, see text. When explicit calibration measurements are 
completed this number could change. 



TABLE IL Parameters in the R-matrix fit of Barker to the beta-delayed a-spectra. 



A 




®B decay 






'^Li decay 








{Mgt)x 


Ex"" 


7A^ 


{Mgt)x 


1 


2.804 


0.588 


0.102 


2.798 


0.591 


0.108 


2 


8.87 


0.884 


-0.180 


8.85 


0.880 


-0.181 


3 


34.1 


1.442 


0.000 


34.1 


1.442 


0.000 





16.72 


0.109 


1.64 


16.72 


0.109 


1.77 


1' 


17.02 


c 


O.OOO'^ 


17.02 


c 


0.000<i 



^This is the channel energy in MeV. 
''In units of MeV^/^^ 

'^a-decay from T = 1 states is isospin forbidden. However, two-state isospin mixing is included, see 
text. 

'^However, a Fermi matrix element of [Mp)ii = \pi is included. 
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FIGURES 

FIG. 1. Recoil corrections, Eqs. ( [B6D and (|B8|), for the electron spectrum of ^Li and the neutrino 
spectrum of ®B respectively, and their ratio. Graph (a) includes the induced pseudotensor current 
form factor; graph (b) excludes the induced pseudotensor current form factor. 



FIG. 2. Radiative corrections, Eqs. and (C9), for the electron spectrum of ®Li and the 

neutrino spectrum of beta decay respectively and their ratio. Graph (a) is the case when the 
internal bremsstrahlung photons are not detected, u large (see Appendix |C|) ; graph (b) is the case 
when the internal bremsstrahlung photons are detected and the total energy deposited, X = E(.+Q, 
is recorded in the SNO detector, to small. 



FIG. 3. Values of the ratio, (Te)^/(Te)e, and their 3a errors due to various uncertainties. Also 
shown are various theoretical values corresponding to different neutrino-oscillation scenarios, taken 
from Bahcall and Lisi ||5|. Labels: STD = standard (no oscillation); SMA = small-mixing angle 
(MSW); LMA = large-mixing angle (MSW); VAC = vacuum oscillation. 

FIG. 4. Ratios of normalized electron spectra as a function of the electron kinetic energy, Tg 
(see text). The two dashed lines enclose the error band without a ®Li calibration experiment, due 
to the ±3(7 uncertainty in absolute energy calibration. The two solid lines enclose the error band 
in the ratio of measured ®B solar-neutrino to measured ^Li spectra under the same assumptions. 
The data points are theoretical values for the CC spectrum evaluated under the small-angle MSW 
solution (SMA). 
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CC — shape uncertainties 
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